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Abstract 



One of the many remarkable features of MHV scattering amplitudes is their conjec- 
tured equality to lightlike polygon Wilson loops, which apparently holds at all orders 
in perturbation theory as well as non-perturbatively. This duality is usually expressed 
in terms of purely four- dimensional quantities obtained by appropriate subtraction of 
the IR and UV divergences from amplitudes and Wilson loops respectively. In this 
paper we demonstrate, by explicit calculation, the completely unanticipated fact that 
the equality continues to hold at two loops through 0{t) in dimensional regulariza- 
tion for both the four-particle amplitude and the (parity-even part of the) five-particle 
amplitude. 



1 Introduction 



Amongst the many remarkable features of the mathematical structure of scattering 
amplitudes that have emerged in the past several years, one of the most mysterious 
remains the apparent equality between planar maximally helicity violating (MHV) 
scattering amplitudes and lightlike Wilson loops in maximally supersymmetric Yang- 
Mills (SYM) theory. This new aspect of duality first emerged in [1], where Alday and 
Maldacena argued using AdS/CFT that the prescription for computing scattering 
amplitudes at strong coupling was mechanically identical to that for computing the 
expectation value of a Wilson loop over the closed contour obtained by gluing the 
momenta of the scattering particles back-to-back to form a polygon with lightlike 
edges. 

Adhering to the principle that there is no such thing as a coincidence in SYM 
theory, it was suggested in [2l[3] that MHV amplitudes and Wilson loops might be 
equal to each other not just at strong coupling, as the work of Alday and Maldacena 
indicated, but perhaps even order by order in perturbation theory. This bold sugges- 
tion was confirmed by explicit calculations at one loop for four particles in [2] and for 
any number of particles in and at two loops for four and five particles in jUE]. 

Already for a few years prior to these developments planar MHV amplitudes in 
SYM had come under close scrutiny following the discovery of the ABDK relation [6] , 
which expresses the four-point two-loop amplitude as a certain quadratic polynomial 
in the corresponding one-loop amplitude, a relation which was later checked to hold 
also for the five-point two-loop amplitude [HIH]. The all-loop generalization of the 
ABDK relation, known as the BDS ansatz after the authors of [9], expresses an 
appropriately defined infrared finite part of the all-loop amplitude in terms of the 
exponential of the one-loop amplitude. This proposal has also been completely verified 
for the three-loop four-point amplitude [H], and partially explored for the three- loop 
five-point amplitude [TD]. 

However it was shown in [llj that the ABDK/BDS ansatz is incompatible with 
strong coupling results in the limit of a very large number of particles, and indeed 
it was found in [12] that starting from six particles and two loops the ansatz is 
incomplete and the amplitude is given by the ABDK/BDS expression plus a nonzero 
'remainder function' (an analytic expression for which was obtained in [T3HI5]). The 
breakdown of the ABDK/BDS ansatz beginning at six particles can be understood 
on the basis of dual conformal symmetry [l6l[T7], which completely determines the 
form of the four- and five-particle amplitudes but allows for an arbitrary function of 
conformal cross-ratios beginning at n = 6 [HITS]. While dual conformal invariance of 
SYM scattering amplitudes remains a conjecture beyond one loop, it is necessary if the 
equality between amplitudes and Wilson loops is to hold in general since the symmetry 
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translates to the manifest ordinary conformal invariance of the corresponding Wilson 
loops. 

Of course dual conformal symmetry alone does not imply the amplitude/ Wilson 
loop equality since they could differ by an arbitrary function of cross-ratios, but 
miraculously precise agreement was found in [T21IIH] between the two sides for n = 6 
particles at two loops. Evidently some magical aspect of SYM theory is at work 
beyond the already remarkable dual conformal symmetry. 

This series of developments has opened up a number of interesting directions 
for further work. In this paper we turn our attention to a question which might 
have seemed unlikely to yield an interesting answer: does the amplitude/Wilson loop 
equality hold beyond 0{e^) in the dimensional regularization parameter e? This 
question is motivated largely by the observation [3] that at one loop, the four-particle 
amplitude is actually equal to the lightlike four-edged Wilson loop to all orders in e 
after absorbing an e-dependent normalization factor. Furthermore, the parity-even 
part of the five-particle amplitude is equal to the corresponding Wilson loop to all 
orders in e, again after absorbing the same normalization factorj^ To our pleasant 
surprise we find a positive answer to this question at two loops: agreement between the 
n = 4 and the parity-even part of the n = 5 amplitude and the corresponding Wilson 
loop continues to hold at 0{e) up to an additive constant which can be absorbed into 
various structure functions. 

Let us emphasize that this is a rather striking result which cannot reasonably be 
called a coincidence: at this order in e the amplitudes and Wilson loops we compute 
depend on all of the kinematic variables in a highly nontrivial way, involving poly- 
logarithmic functions of degree 5. It would be very interesting to continue exploring 
this miraculous agreement and to understand the reason behind it. Dual conformal 
invariance cannot help in this regard since the symmetry is explicitly broken in di- 
mensional regularization so it cannot say anything about terms of higher order in e, 
but of course as mentioned above already at 0{e^) there must be some mechanism 
beyond dual conformal invariance at work. 

The rest of the paper is organized as follows. In Section 2 we review key aspects 
of one- and two- loop amplitudes and Wilson loops, in particular the ABDK/BDS 
ansatz and the correspondence between MHV amplitudes and Wilson loops. We also 
summarize our main results on the equality of amplitudes and Wilson loops up to 0(e) 
at four and five points in fl2.17p . f l2.18p . respectively. In Section 3 we present the four- 
and five-point amplitudes at one and two loops, and in Section 4 the corresponding 
Wilson loops. Section 5 is devoted to discussing the numerical methods that have 

-'^For n > 5 the Wilson loop calculation reproduces only the all orders in e two-mass easy box 
functions, while the corresponding n-point amplitude contains additional parity-odd as well as parity- 
even terms which vanish as e — 0. 
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been employed in order to perform our analysis. Finally, in Section 6 we compare 
amplitude and Wilson loops, showing the agreement between these two quantities up 
to and including C(e) terms. Two appendices complete the paper. In the first one, 
we present results valid to all orders in e for all Wilson loop diagrams contributing 
to the four-point case, with the exception of the so-called "hard" diagram topology, 
which is evaluated only up to and including C(e) terms. In the second appendix, 
we present a novel expression for the all-orders in e one-loop ra-point Wilson loop 
diagrams which have simple analytic continuation properties. 



2 Review and Summary of Main Results 



The infinite sequence of n-point planar maximally helicity violating (MHV) ampli- 
tudes in A/" = 4 super- Yang-Mills theory (SYM) has a remarkably simple structure. 
Due to supersymmetric Ward identities [T91422] . at any loop order L, the amplitude 
can be expressed as the tree-level amplitude, times a scalar, helicity-blind function 

4!^^ =ATMl!^l (2.1) 

In [B], ABDK discovered an intriguing iterative structure in the two- loop expansion 
of the MHV amplitudes at four points. This relation can be written as 

M?{e)-^-{M^\e)Y = /(^)(e)A^W(2e) + C(^) + 0(e) , (2.2) 

where IR divergences are regulated by working in D = 4 — 2e dimensions (with e < 0), 

f^'\e) = -C-Cse-Ue' , (2.3) 

and 

C(^) = -^Cl . (2.4) 

The ABDK relation (12. 2 p is built upon the known exponentiation of infrared diver- 
gences [231124] , which guarantees that the singular terms must agree on both sides of 
(12. 2p . as well as on the known behavior of amplitudes under collinear limits |25[|26] . 
The (highly nontrivial) content of the ABDK relation is that (12. 2 p holds exactly as 
written at 0(e°). However, ABDK observed that the 0{e) terms do not satisfy the 
same iteration relation 0. 

In [6], it was further conjectured that (12. 2p should hold for two-loop amplitudes 
with an arbitrary number of legs, with the same quantities (12. 3p and (12. 4p for any n. 
In the five-point case, this conjecture was confirmed first in [7] for the parity-even part 
of the two-loop amplitude, and later in [8J for the complete amplitude. Notice that 
for the iteration to be satisfied parity-odd terms that enter on the left-hand side of 
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the relation must cancel up to and including C(e°) terms, since the right-hand side is 
parity even up this order in e. So far this has been checked and confirmed at two-loop 
order for five and six particles [8l|27]. This is also crucial for the duality with Wilson 
loops (discussed below) which by construction cannot produce parity-odd terms at 
two loops. 



It has been found that starting from six particles and two loops, the ABDK/BDS 
ansatz (12.21) needs to be modified by allowing the presence of a remainder function 

7^n P/18J, 

M^:\e) - i(A^W(e))' = f^'\e)M^^\2e) + C^'^ + 7^„ + E^{e) , (2.5) 
where 7^„ is e- independent and En vanishes as e — )■ 0. We parameterize the latter as 

En{e) = tSn + 0{t^) . (2.6) 

In this paper we will discuss in detail for n = 4, 5 where we find a remarkable 
relation to the same quantity calculated from the Wilson loop. Hitherto this relation 
was only expected to hold for the finite parts of the remainder 7^„. 



In a parallel development, Alday and Maldacena addressed the problem of calcu- 
lating scattering amplitudes at strong coupling in A/" = 4 SYM using the AdS/CFT 
correspondence. Their remarkable result showed that the planar amplitude at strong 
coupling is calculated by a Wilson loop 



W[Cn] := TrPexp 



(2.7) 



whose contour C„ is the ?7,-edged polygon obtained by joining the lightlike momenta 
of the particles following the order induced by the colour structure of the planar 
amplitude. At strong coupling the calculation amounts to finding the minimal area 
of a surface ending on the contour C„ embedded at the boundary of a T-dual AdS^ 
space [1]. Shortly after, it was realised that the very same Wilson loop evaluated at 
weak coupling reproduces all one-loop MHV amplitudes in A/" = 4 SYM [2l|3]. The 
conjectured relation between MHV amplitudes and Wilson loops found further strong 
support by explicit two loop calculations at four , five [5] and six points [T^IIISIEZ] • 
In particular, the absence of a non-trivial remainder function in the four- and five- 
point case was later explained in ^ from the Wilson loop perspective, where it was 
realised that the BDS ansatz is a solution to the anomalous Ward identity for the 
Wilson loop associated to the dual conformal symmetry [16]. 



The Wilson loop in (12. 7p can be expanded in powers of the 't Hooft coupling 
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a := g^N/{8TT^) at| 



{W[Cn]) := 1 + J^a'iyW := exp^a'^/;« . (2.8) 
1=1 1=1 

Note that the exponentiated form of the Wilson loop is guaranteed by the non-Abelian 
exponentiation theorem [29||30] . The Wn^ are obtained from "maximally non-Abelian" 
subsets of Feynman diagrams contributing to the Wn'^ and in particular from fl2.8p 
we find 

n^':' = , ^1^^ = Wi^^ - \ (W^W)^ . (2.9) 

The UV divergences of the ra-gon Wilson loop are regulated by working in D = 4 + 2e 
dimensions with e < 0. The one-loop Wilson loop Wn^ times the tree-level MHV 
amplitude is equal to the one-loop MHV amplitude, first calculated in [31] using the 
unitarity-based approach [32], up to a regularization-dependent factor. This implies 
that non-trivial remainder functions can only appear at two and higher loops. At two 
loops, which is the main focus of this paper, we define the remainder function Tl^^ 
for an n-sided Wilson loop a^ 



<K^) = P4l{^) KK'^^) + c'4i + T^T + , (2.10) 



where 



:= /o^'^ + ffk' + f?k^ ■ (2.11) 
Note that /q = —(^2, which is the same as on the amplitude side, while ~ 

/ox 

^eik — '^Cs [33j. In t28j, the four- and five-edged Wilson loops were cast in the form 
f l2.10p and by making the natural requirements 

^WL ^ ^WL ^ Q ^ ^2.12) 

this allowed for a determination of the coefficients /2 ^wl- The results found 
in [2H], ar(EI 

fSii^) = -C2 + 7Cse - 5C,e' , (2.13) 

and 

C^i = . (2.14) 

As noticed in ^28j, there is an intriguing agreement between the constant ^^'1 
the corresponding value of the same quantity on the amplitude side. 



^We follow the definitions and conventions of [28|, to which we refer the reader for more details. 
■^We expect a remainder function at every loop order I and the corresponding equations would 



*The 0(1) and 0(e) coefficients of /^lI^) h^'d been determined earlier in [4]. 
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What has been observed so far is a duahty between Wilson loops and amplitudes 
up to finite terms. In turn this can be reinterpreted as an equality of the corresponding 
remainder functions [f| 

7^„ = 7^^ • (2.15) 

(2) (2) 

A consequence of the precise determination of the constants /2 WL ^"^^ ^wl ^^^^ 
no additional constant term is allowed on the right hand side of fl2.15p . For the same 
reason, the Wilson loop remainder function must then have the same coUinear limits 
as its amplitude counterpart, i.e. 

7^^L ^ n'^^ , (2.16) 

with no extra constant on the right hand side of (I2.16p . 

The main result of the present paper is that for n = 4, 5 the relation between 
amplitudes and Wilson loops continues to hold for terms of order e^. In particular we 
find 

= ^i'wL - 3C5 , (2.17) 

4"^=4%L-IC^- (2.18) 

Note that these results have been obtained (semi-)numerically with typical errors of 
10~^ at n = 4 and 10~^ for n = 5. Details of the calculations are presented in the 

remaining sections of this paper. More precisely S^^ is known analytically [H], while 

(2) 

the analytic evaluation of discussed in appendix A. At five points all results 

are numerical and furthermore on the amplitude side we only considered the parity- 
even terms. It is an interesting open question whether the parity-odd terms cancel 
at C(e) as they do at C(e°) [8]. 



3 Amplitudes 

In this section we review the ingredients necessary for our calculation of the 0{e) 
terms in the ABDK relation for the n = 4, 5 point amplitudes. 



3.1 One-Loop Amplitudes 

We begin with the one-loop amplitudes, for which analytic results can be given to all 
orders in e. 

^An alternative interpretation of the duality in terms of certain ratios of amplitudes (Wilson 
loops) has been given recently in f M] . 
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Following the conventions of [6J, the one- loop four-point amplitude may be ex- 
pressed as [35] 



M 



(1) 



— St I 



(1) 



(3.1) 



where s = {pi + ^2)^, t = {p2 + Ps^ are the usual Mandelstam variables and I^^^'' is 
the massless scalar box integral 

P2 P3 



r(l) 



3^7 



ITT 



D/2 



d^p 



p^{p - piY{p -pi- P2Y{p + p^y 



(3.2) 



Pi 



Pa 



which we have written out in order to emphasize the normalization convention (fol- 
lowed throughout this section) that each loop momentum integral carries an overall 
factor of e'^'^/ivr^/^. The integral may be evaluated explicitly (see for example [36]) in 
terms of the ordinary hypergeometric function 2F1, leading to the exact expression 



M 



(1) 



e^^r(l + e)r^(l-e) 
r(l-2e) 



[(-s)-^2Fi(l, -e, 1 - e, 1 + s/t) + {s^ t)] , (3.3) 



valid to all orders in e. We will always be studying the amplitude/ Wilson loop 
duality in the fully Euclidean regime where all momentum invariants such as s and 
t are negative. The formula (13. 3p applies in this regime as long as we are careful to 
navigate branch cuts according to the rule 



-2)~'2Fi(-e, -e, 1 - e, 1 + z) := limRe 

e->-0 



2Fi(-e, -e,l - e,l + z + is) 
{-z + ieY 



(3.4) 



when z > 0. 



Five-point loop amplitudes A^s^"* contain both parity-even and parity-odd contri- 
butions after dividing by the tree amplitude as in (12. ip . The parity-even part of the 
one-loop five-point amplitude is given by [3j 



P3 



Pi 



M 



(1) 

5+ 



(1) 
5 ' 



r(l) 



cyclic 



P2 



Pi 



Pb 



(3.5) 



where Sj = {pi+Pi+i)"^ and the sum runs over the five cyclic permutations of the exter- 
nal momenta Pi. This integral can also be explicitly evaluated (see for example [36]). 
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Figure 1: Integrals appearing in the amplitude A^5+- Note that I^^^'^ contains the 
indicated scalar numerator factor involving q, one of the loop momenta. 

leading to the all-orders in e result 

r(l - 2e) 2 ^ l\ S3S4 / si - S4 

cyclic ^ ' 

■Si — S3 V T7 / 1 1 "^4 

2F1 -e, -e, 1 -e, 1 



S3S4 / Si- S3 

(Si - S3)(si - S4) V „ . ^ ^ S3S4 

2ti[-e, -e, i - e, i 



(3.6) 

again keeping in mind (I3l 



3.2 Two-Loop Amplitudes 



The two-loop four-point amplitude is expressed as [38] 

P2 



P3 



A^f) = ]sHlP + {s^t), 



r(2) 



(3.7) 



Pi 



Pi 



which may be evaluated analytically through 0{e^) using results from [9j (no all-orders 
in e expression for the double box integral is known), from which we find 



^4 = 5 Li5(-x) - 4LLi4(-a;) + ^(3^^ + tt^) Li3(-a;) - ^(L^ + vr^) Li2(-x) 
- 1(L^ + n^f log(l + a:) + l^rU - ^(5 + gvr^Cs, 



(3.^ 



where x = t/s and L = logx. A comment is in order here: In order to be able to 
present the amphtude remainder (13.81) in this form, we have pulled out a factor of 
{st)~^^^'^ from each loop amplitude This renders the amplitudes, and hence the 

ABDK remainder E^le), dimensionless functions of the single variable x. We perform 
this step in the four-point case only, where we are able to present analytic results for 
the amplitude and Wilson loop remainders. 

The parity-even part of the two-loop five-point amplitude involves the two integrals 
shown in Figure [H in terms of which [Tf [81139] 

= ^ E (^3^^/^'^" + {P^ ^ Pe-^)) + S.SsSj^'^', (3.9) 
cyclic 

where Si = {pi + Pi+i)^. To evaluate this amplitude to C(e) we must resort to a 
numerical calculation using Mellin-Barnes parameterizations of the integrals (which 
may be found for example in [^), which we then expand through 0{e), simplify, 
and numerically integrate with the help of the MB, MBresolve, and barnesroutines 
programs pO]HT| . In this manner we have determined the 0{e) contribution S^^^ to the 
five-point ABDK relation numerically at a variety of kinematic points. The results 
are displayed in Table [H 



4 Wilson Loops 



4.1 One-Loop Wilson Loops 

The one-loop Wilson loop was found in [3] for any number of edges and to all orders in 
the dimensional regularization parameter e. It is obtained by summing over diagrams 
with a single gluon propagator stretching between any two edges of the Wilson loop 
polygon. Diagrams with the propagator stretching between adjacent edges pi and 




cusp :=r(l + e)e'^ x {-^{si)-') finite :=r(l + e)e'^ x J", 



Figure 2: One-loop Wilson loop diagrams. The expression of is given in ( IB.12P of 
equivalently in (lB.14p . 

Pi+i are known as cusp diagrams, and give the infrared- divergent terms in the Wilson 
loop, proportional to {—2pi ■ pi^i)~'' / = (— Sj)~Y^^- 
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On the other hand, diagrams for which the propagator stretches between two non- 
adjacent edges are finite. Their contribution to the Wilson loop can be found to all 
orders in e and is (up to an e-dependent factor) precisely equal to the finite part of a 
two- mass easy or one-mass box function [3] (for details see appendix [B]) . The general 
rz-point one loop amplitude is given by the sum over precisely these two-mass easy 
and one-mass box functions [S] to 0{e^)E Thus we conclude that the Wilson loop 
is equal to the amplitude at one loop for any n up to finite order in e only (and up 
to a kinematic independent factor). 



However at four and five points a much stronger statement can be made. The 
four-point amplitude and the parity-even part of the five-point amplitude are both 
given by the sum over zero- and one-mass boxes to all orders in e. Thus the Wilson 
loop correctly reproduces these one-loop amplitudes to all orders in e. Using the 
results in appendix B we find that the four-point Wilson loop (in a form which is 
manifestly real in the Euclidean regime s,t < 0) is given by 



W, 



(1) 



r(l + e)e^^| - ^ [i-s)-^ + i-t)-'] + T,{s, t, 0, 0) + J^,{t, s, 0, 0)} 



e2 \st. 



2Fi(e,e;l + e;-t/s)+ (^^^ 2Fi(e, e; 1 + e; -s/t) - 27re cot(e7r) 

(4.1) 



Note in particular the additional cotangent term explained in detail at the end of 
appendix B. The generic form of the function J-'^ is given in flB.12p of equivalently in 

dEHD. 



For the five-point amplitude we display a new form which has a simple analytic 
continuation in all kinematical regimes and also a very simple expansion in terms 
of Nielsen polylogarithms (see (IB. lip ). It is given in terms of 3F2 hypergeometric 



^The all-orders in e n-point amplitude contains new integrals contributing at 0(e). 
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functions and is derived in detail in appendix B: 

5 



Er(i + 



e e 



i=l 



+ 



1 / -Sj+S ^ ■5j — Sj-|_i 



5?3F2(l,l,l + e;2,2;^) 
+ ^'^^ 3F2 (1, 1, 1 + e; 2, 2; 



(4.2) 



where Hn is the n^'^-harmonic number. Using hypergeometric identities one can show 
that (up to the prefactor) the four- and five-sided Wilson loops (14. ip . (14. 2 p are equal 
to the four-point and the (parity-even part of the) five-point amplitudes of (13. 3p 
and (EH). 



The precise relation between the Wilson loop and the amplitude is 



(1) _ r(i-26) 



(1) 



where A^4^'* is the one-loop four-point amplitude and A^s!/ is the parity-even part of 
the five-point amplitude. 



r(i-26) 



(1) 



(4.3) 



4.2 Two-Loop Wilson Loops 



At two-loop order, the n-point Wilson loop is given by a sum over six different types 
of diagrams. These are described in general for polygons with any number of edges 
in [28] and are displayed for illustration below. 

The computation of the four-point two-loop Wilson loop up to C(e°) was first 
performed in [1]. In appendix |X] we display all the contributing diagrams for this 
case and give expressions for these to all orders in e in all cases except for the "hard" 
diagram, which we give up to and including terms of 0{e). Summing up the con- 
tributions from all these diagrams we obtain the result for the two-loop four-point 
Wilson loop to 0{e). This is displayed in (14. 4 p of the next subsection. 
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cross curtain factorised cross 




Y hard self-energy 



Figure 3: The six different diagram topologies contributing to the two-loop Wilson 
loop. For details see f2Bj . 

The five-point two-loop Wilson loop was calculated up to 0(e°) in [5]. In order to 
obtain results at one order higher in e we have proceeded by using numerical methods. 
In particular we have used Mellin-Barnes techniques to evaluate and expand all the 
two-loop integrals of Figure [3l This is described in more detail in Section [5l 



4.2.1 The Complete Two-Loop Wilson Loop at Four Points 



Here is our final result for the four-point Wilson loop at two loops expanded 
and including terms of 0{e): 



w 



(2) 



C X 



{-s)-'' + {-t) 



-2t 



where 

W2 = 

Wi = 

Wo = 

W-i = 



48 



rC3 



-— (log^a; + TT^) + 
48 ^ ^ ^144 



TT 



log X H TT 

48 V ^ 3 



up to 
(4.4) 

(4.5) 
(4.6) 
(4.7) 



1440 



467r'^logx- lOvr^ \og^ x + 75tc'^ log(l + x) + 90%^ \og^ x \og{l + x) 



-M51og^xlog(l +x) + 2407r^logxLi2(-x) + 120 log^ x Li2(-x) 
-30077^113 (-a;) - 540 log^ x Li3(-x) + 1440 logx Li4(-a;) 

-1800Li5(-x) - 156077^(3 - 12601og2a;C3 + 5940(5 



(4. 
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and 



C := 2 



r(l + e)e 



7e 



2(l + C26^-^C3e^) +0{e' 



(4.9) 



We recall that x = t/s. 



We would like to point out the simplicity of our result fl4.4p - specifically, f l4.5p - 
(14. 8p are expressed only in terms of standard polylogarithms. Harmonic polyloga- 
rithms and Nielsen polylogarithms are present in the expressions of separate Wilson 
loop diagrams, as can be seen in appendix A, but cancel after summing all contribu- 
tions. 



4.2.2 The 0(e) Wilson Loop Remainder Function at Four Points 



Using the result ( 14. 4p and the one-loop expression for the Wilson loop, one can work 
out the expression for the remainder function at 0{e), as defined in (12. 5p and (12. 6p . 
Our result is 



4,WL 



1 

360 



IQir"^ log X — IStt"^ log(l + x) — 307r^ log^ x log(l + x) 



-151og^xlog(l + x) - 1207rMogxLi2(-x) - 120 log^(x)Li2(-x) 



+1807r2Li3(-x) + 540 log^ xLi3(-x) 
+1800Li5(-x) + 69071^3 - 5940(5 



14401og(x)Li4(-a;) 



(4.10) 



where we recall that £n,WL is related to the quantity En introduced in (12.50 and (12. 6p . 
Remarkably, (I4.10p does not contain any harmonic polylogarithms. We will compare 
the Wilson loop remainder (14.1 OP to the corresponding amplitude remainder (13.80 in 
Section Elfl 



4.2.3 The (9(e) Wilson Loop at Five Points and the Five-Point Remainder 
Function 



For the five-point amplitude and Wilson loop at two loops we resort to completely 
numerical evaluation of the contributing integrals, and a comparison of the remainder 
functions is then performed. We postpone this discussion to section 16.21 

^ Similarly to what was done for the amplitude remainder p.Sp . in arriving at (|4.10|) we have 
pulled out a factor of (st)^*^/^ per loop in order to obtain a result which depends only on the ratio 
X := t/s. 
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5 Mellin-Barnes Integration 



The two-loop five-point Wilson loop and amplitude have been numerically evalu- 
ated by means of the Mellin-Barnes (MB) method using the MB package [40] in 
MATHEMATICA. At the heart of the method lies the Mellin-Barnes representation 



/+ioo -y-z 
dz-^Ti-znX + z). (5.1) 
■too 



(X + Y)^ 2TTiT{X) 

We will use the integral representation for the hard diagram of the Wilson loop as an 
example in order to describe the procedure we followed. The integral for the specific 
diagram shown in Figure H] has the expression 

/h(Pi,P2,P3;Qi,Q2,Q3) (5.2) 
_ ir(2 + 2e) ^e AT 



rri+ )2 / (Yl^^i) / (Yldai)5{l-^ai){aia2a3) 



-2e ' 



We write the numerator and denominator as a function of the momentum invariants, 
i.e. squares of sums of consecutive momenta, 

V = -a^a2 [{pi + Qs + P2?{1 - ri)T2 + {pi + Q^Yil - n){l - T2) 

+ {Q3+P2fnT2 + Qlnil-T2)]+cyc\ic{l,2,3), (5.3) 



J\f =2 [2{piP2)ip3Q3) - iP2P3)iPlQ3) - {PlP3){P2Q3)] 011^2 

+ 2{piP2){pzPi) [ai«2(l - Ti) + ttsaiTi] + cyclic(l, 2, 3), (5.4) 

where 

2piPi+i = -{Pi + Qi+2f + Ql+2 - {Qi+2 + Pi+if + {Qi + Pi+2 + Qi+if, 
'^PiQi = -{pi + Qi+2 + Pi+if + {Qi+2 + Pi+if 

- {pi+2 + Qi+1 + Pif + {Pi+2 + Qi+lf, 
'2p^Q, = {p^ + Q,f -Q]. (5.5) 

By means of the substitution ai — )■ 1 — r4, a2 — )• T/^t^ and — )■ r4(l — ts), we eliminate 
one integration and the delta function to get a five-fold integral over Xj G [0, 1]. Next, 
we obtain an MB representation using the generalisation of fl5.ll) 




(e:LiX.)a {2mr-^T{\) 



which introduces m— 1 MB integration variables Zg-, where m is the number of terms in 
the denominator. At this point, the integrations over the r^'s can be easily performed 
by means of the substitution 
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Pi 



Pi 




P2 



Figure 4: The hard diagram corresponding to fl5.2p . 



We are now left with an integrand that is an analytic function containing powers 
of the momentum invariants {—Sij)'^^^^''^'^^ and Gamma functions T(g({zi}, e)), where 
/ and g are linear combinations of the z^^s and e. In order to perform the MB 
integrations, one has to pick appropriate contours, so that for each Zg the r(- ■ ■ + Zg) 
poles are to the left of the contour and the r(- ■ ■ — z^) poles are to the right. 



At this point we use various Mathematica packages to perform a series of op- 
erations in an automated way to finally obtain a numerical expression at specific 
kinematic points. We will briefly summarise the steps followed, while for more de- 
tails we refer the reader to the references documenting these packages and references 
therein. Using the MBresolve package [51], we pick appropriate contours and resolve 
the singularity structure of the integrand in e. The latter involves taking residues 
and shifting contours, and is essential in order to be able to Laurent expand the inte- 
grand in e. Using the barnesroutines package [^[¥T]. we apply the Barnes lemmas, 
which in general generate more integrals but decrease their dimensionality, leading to 
higher precision results. Finally, using the MB package PO] we numerically integrate at 
specific Euclidean kinematic points to obtain a numerical expression. While all ma- 
nipulations of the integrals and the expansion in e are performed in Mathematica, the 
actual numerical integration for each term is performed using the CUBA routines [12] 
for multidimensional numerical integration in FORTRAN. The high number of diagrams, 
and number of integrals for each diagram, makes the task of running the FORTRAN 
integrations ideal for parallel computing. 
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6 Comparison of the Remainder Functions 



6.1 Four-point Amplitude and Wilson Loop Remainders 

The remainder functions for the four-point amphtude and Wilson loops are given in 
(I3.8p and fl4.10p . respectively. From these relations, it follows that the difference of 
remainders is a constant, x- independent term: 



as anticipated in f l2.17p . 

We would like to stress that this is a highly nontrivial result since there is no 
reason a priori to expect that the four-point remainder on the amplitude and Wilson 
loop side, (13 -Sp and fl4.10p respectively, agree (up to a constant shift). For example, 
anomalous dual conformal invariance is known to determine the form of the four- and 
five-point Wilson loop only up to 0{e^) terms [5], but does not constrain terms which 
vanish as e — )■ 0. The expressions we derived for the amplitude and Wilson loop 
four-point remainders at 0{e) are also pleasingly simple, in that they only contain 
standard polylogarithms. 



6.2 Five-point Amplitude and Wilson Loop Remainders 

We have numerically evaluated both the five-point two-loop amplitude and Wilson 
loop up to C(e) at 25 Euclidean kinematic points, i.e. points in the subspace of the 
kinematic invariants with all Sij < 0. The choice of these points and the values 
of the remainder functions £^ , ^^s^l ^(^) together with the errors reported by 
the CUBA numerical integration library [12] appear in Table [H while in Figures 
|5] and |6] we plot both remainders for all kinematic points. We have calculated the 
difference between the amplitude and Wilson loop remainders, see Table [2] and Figure 
[71 Remarkably, this difference also appears to be constant (within our numerical 
precision) as in the four-point case, and hence we conjecture that 



It is also intriguing that the constant difference is fit very well by a simple rational 
multiple of Cs, rather than a linear combination of and C2C3 as would have been 
allowed more generally by transcendentality. 

In the last column of Table[2]we give the distance of our results from this conjecture 
in units of their standard deviation. 



S4 — £'4,WL — 3^5 , 



(6.1) 



■(2) _ p(2) 
5 ~ ^5,WL 




(6.2) 
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# 


(Sl2, S23, 


■534, 


945, ssi) 




<r(2) 
5,WL 


1 


(-1,-1, 


-1, 


-1,-1) 


-8.463173 ± 0.000047 


-5.8705280 ± 0.0000068 


2 


(-1,-1, 


-2, 


-1,-1) 


-8.2350 ± 0.0024 


-5.64560 ± 0.00063 


3 


(-1,-2, 


-2, 


-1,-1) 


-7.7697 ±0.0026 


-5.17647 ±0.00076 


4 


(-1,-2, 


-3, 


-4, -5) 


-6.234809 ±0.000032 


-3.642125 ± 0.000018 


5 


(-1,-1, 


-3, 


-1,-1) 


-8.2525 ±0.0027 


-5.65919 ±0.00097 


6 


(-1,-2, 


-1, 


-2,-1) 


-8.142702 ±0.000023 


-5.5500050 ± 0.0000092 


7 


(-1,-3, 


-3, 


-1,-1) 


-7.6677 ±0.0034 


-5.0784 ±0.0013 


8 


(-1,-2, 


-3, 


-2,-1) 


-6.8995 ± 0.0029 


-4.31395 ± 0.00093 


9 


(-1,-3, 


-2, 


-5, -4) 


-6.9977 ±0.0031 


-4.40806 ± 0.00099 


10 


(-1,-3, 


-1, 


-3,-1) 


-8.2759 ±0.0025 


-5.69086 ± 0.00085 


11 


(-1,-4,- 


-8, - 


16, -32) 


-8.7745 ±0.0078 


-6.1825 ±0.0051 


12 


(-1,-8,- 


-4, - 


32, -16) 


-11.991985 ±0.000089 


-9.398659 ± 0.000084 


13 


(-1, -10, - 


-100 


-10,-1) 


-2.914 ±0.022 


-0.300 ±0.010 


14 


(-1, -100, - 


-10, 


-100, -1) 


-3.237 ±0.011 


-0.6648 ± 0.0028 


15 


(-1,-1,- 


-100 


-1,-1) 


-12.686 ± 0.014 


-10.108 ± 0.010 




(A 1 on 


-1, 


-100, -1) 


1 A 7f\R7 -t- n 0077 
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17 


(-1, -100, 


-100,-1,-1) 


-182.32 ±0.11 


-179.722 ± 0.039 


18 


(-1, -100, - 


-10, 


-100, -10) 


-6.3102 ±0.0062 


-3.7281 ±0.0013 


19 
20 
21 


(-1,-1, 


1 

r 

*1 
1' 

4 ' 


— -— ) 

¥' ¥ 

25' 16/ 
-1,-1) 


-19.0031 ±0.0077 
-15.1839 ±0.0046 
-9.7628 ±0.0028 


-16.4136 ±0.0021 
-12.5995 ±0.0016 
-7.17588 ±0.00079 


22 


(-i.-f 
(-i.-i 


-1,-1) 


-9.5072 ± 0.0036 


-6.9186 ±0.0014 


23 


-1, 


-^'-1) 


-12.6308 ±0.0031 


-10.04241 ±0.00083 


24 


(-1,-i 


1 

?' 

4' 


-11.0200 ±0.0056 


-8.4281 ± 0.0030 


25 


(-1,-i 




-19.1966 ±0.0070 


-16.6095 ±0.0043 



Table 1: C(e) five-point remainders for amplitudes (£^5 ) and Wilson loops (£^5 wl)- 



# 


(si2, S23, 


S34, S45, S51) 


•--5 '-•5,WL 




1 


(-1,-1, 


-1,-1,-1) 


-2.592645 ± 0.000048 


6.8 


2 


(-1,-1, 


-2,-1,-1) 


-2.5894 ±0.0025 


1.2 


3 


(-1,-2, 


-2,-1,-1) 


-2.5932 ±0.0027 


0.32 


4 


(-1,-2, 


-3,-4,-5) 


-2.592697 ± 0.000036 


10 


5 


(-1,-1, 


-3,-1,-1) 


-2.5933 ± 0.0028 


0.35 


6 


(-1,-2, 


-1,-2,-1) 


-2.592697 ± 0.000025 


15 


7 


(-1,-3, 


-3,-1,-1) 


-2.5893 ± 0.0036 


0.82 


8 


(-1,-2, 


-3,-2,-1) 


-2.5856 ± 0.0030 


2.2 


9 


(-1,-3, 


-2,-5, -4) 


-2.5897 ±0.0032 


0.82 


10 


(-1,-3, 


-1,-3,-1) 


-2.5851 ±0.0026 


2.8 


11 


(-1,-4,- 


-8, -16, -32) 


-2.5920 ± 0.0093 


0.034 


12 


(-1,-8,- 


-4, -32, -16) 


-2.59333 ±0.00012 


8.3 


13 


(-1, -10, - 


-100,-10,-1) 


-2.614 ±0.024 


0.89 


14 


(-1, -100, - 


-10, -100, -1) 


-2.572 ±0.011 


1.9 


15 


(-1,-1,- 


-100, -1, -1) 


-2.578 ±0.017 


0.80 


16 


(-1, -100, 


-1,-100,-1) 


-2.593 ±0.010 


0.071 


17 


(-1, -100, 


-100,-1,-1) 


-2.60 ±0.11 


0.039 


18 


(-1, -100, - 


-10, -100, -10) 


-2.5820 ± 0.0063 


1.6 


19 


(-1.-I.- 


- -— - — ) 

?' 2^5 
4 ' 25 ' 16) 


-2.5894 ±0.0080 


0.36 


20 


-2.5844 ± 0.0049 


1.6 


21 


(-1,-1, 


4.-1.-1) 


-2.5869 ±0.0029 


1.9 


22 


(-l.-i 

(-I'-i 


-2.5886 ±0.0038 


0.96 


23 


-l.-i.-l) 
^ ^ 1) 

-i'-i'- 

4 . q ' ^) 


-2.5884 ±0.0032 


1.2 


24 


(-1,-1, 


-2.5919 ±0.0064 


0.064 


25 


-2.5870 ±0.0082 


0.65 



Table 2: Difference of the five-point amplitude and Wilson loop two-loop remainder 
functions at 0{e), and its distance from — |C5 ~ —2.592319 in units of a, the standard 
deviation reported by the CUBA numerical integration package l4^ . 
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Figure 5: Remainder functions at 0{e) for the amplitude (circle) and the Wilson loop 
(square). 
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Kinematic point 
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Figure 6: Remainder functions at 0{e) for the amplitude (circle) and the Wilson 
loop (square). In this Figure we have eliminated data point 17 and zoomed in on the 
others. 
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For kinematic points 1, 4 and 6 we have evaluated the remainder functions with 
even higher precision, and found agreement with the conjecture to 4 digits. A remark 
is in order here. By increasing the precision, the mean value of the difference of 
remainders approaches the conjectured value, but we notice that in units of a it 
drifts away from it, hinting at a potential underestimate of the errors. To test our 
error estimates we used the remainder functions at 0(e°), that are known to vanish. 
Our analysis confirmed that, as we increase the desired precision, the actual precision 
of the mean value does increase, but on the other hand reported errors tend to become 
increasingly underestimated. 

-2.40r 



-2.45 
-2.50 



to 
I 

in 

to 



-2.55 
-2.60 
-2.65 



_i 1. 



-J-i- 



-2.70 



-2.75' 

5 10 15 20 25 

Kinematic point 

Fi gure 7: Difference of the remainder functions si^ ^ — ^l- 
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A Details of the Two-Loop Four- Point Wilson Loop 
to All Orders in e 



In this appendix, we present the resuhs for the separate classes of Wilson loop di- 
agrams contributing to a four-point loop. In all cases (with the exception of the 
"hard" diagram) our results are valid to all orders in the dimensional regularization 
parameter e. These expressions are given in terms of hypergeometric functions. We 
also expand these up to 0{e) with the help of the mathematica packages HPL and 
HypExp pilil]. 



A.l Two-loop Cusp Diagrams 



Figure 8: The two-loop cusp corrections. The second diagram appears with its mirror 
image where two of the gluon legs of the three-point vertex are attached to the other 
edge; these two diagrams are equal. The blue bubble in the third diagram represents 
the gluon self- energy correction calculated in dimensional reduction. 

The total contributions of all diagrams that cross a single cusp (the crossed dia- 
gram across a cusp, the self-energy diagram across the cusp and the vertex across the 
cusp, as depicted in Figure [H]) is easily seen to b^ 

'r(l + 2e)r(l-e) 



\-2e 



16e4 



r(i + 6) 



Adding the contributions for the four cusps, we obtain 



cusp 



-2e 



i-t) 



-2t 



r(l + 2e)r(l -e) 



Le2 24 



-2e 



(A.l) 



(A.2) 



(A.3) 



*In this and the following formulae, a faetor of C is suppressed in each diagram, where C is defined 



20 



Pi 




Figure 9: One of the four curtain diagrams. The remaining three are obtained by 
cyclic permutations of the momenta. 

A. 2 The Curtain Diagram 



The contribution of all four curtain diagrams is 



curtain 



ni 



\-2e 



-t)-^' 

e(-^C3log^x + ^C5-gc3)+0(e^) 



r(l-2e) 

2 



Ik 

e 2 



71 

160 



TT 

48 



log^ X 



(A.4) 



A. 3 The Factorised Cross Diagram 



The factorised cross diagram is given by the product of two finite one-loop Wilson 
loop diagrams, expressed each by (14. ip 



1 fst 
2^ W 



X ^2Fi(e, e, 1 + e, -l/x) + a;^2Fi(e, e, 1 + e, -x) - 27recot(e7r) 



The result for the factorised cross is therefore 

-2e 



1 /.St 



Fi(e, e; 1 + e; — + 2F1 ( e, e; 1 + e; ] x ^ — 27re cot(e7r 



X 



1 2 



(A.5) 
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P3 



P2 




Pi 



Figure 10: One of the four factorised cross diagram. 



^2 ki k2 h 

Pi Pi 




Pi 



Pi P2 



P3 



Pi 



Figure 11: The Y diagram together with the self-energy diagram. The sum of these 
two topologies gives a maximally transcendental contribution. 



with 



9o 
9-1 



-^(log^x + vr^)' , 



log X + tt' 



log'^a; — 6 log(a; + 1) log^ x — 12Li2(— x) logx 



192 

+ Svr^loga; - Gtt^ log(x + 1) + 12 Lis (-x) - 12(3 



(A.6) 



(A.7) 
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A. 4 The Y Diagram 



The diagrams in Figure [TT] correspond to the following contribution to the two-loop 
Wilson loop: 



1 r(l + 2e) 



16e r(l + e)2 
•1 pi 

X / dTi / dT2 

'0 Jo 



u 



(A. 



2/(zi(ri),Z2(T-2), 2:2(^2)) - I{,Zi{ri),Z2{T2),ki) - I{Zi{ri),Z2{T2),k2) 

where ^i(ri) = /C3 -psn, Z2{t2) = h - P1T2, and 



I{Zi,Z2,Z3) 



da{a{l - a)y - {-Zi + az2 + (1 - (t)2;3)^ + «0 



(A.9) 



where 2^2,^3 niust be lightlike. The evaluation of ( 1A.8I) gives 



st\-^' 1 



+ 4x 



u J 64e4 
r(-e + l)^ 
r(-2e+l) 



2(x+ 1) 



_2,r(l + 2e)r(-e+l) 



r(i + 6) 

2Fi(e, 1 + 2e; 1 + e; -a;) + x^' 2^1 (2e, 2e; 1 + 2e; -a;) 



- 47re cot(27re) + r(l + 2e)r(-2e + 1) + - 

We multiply by four to obtain the contribution of all such diagrams. Then the 
expansion of this contribution in e begins at (9(e"^), 



+ i-ty 



X 



Cl 



+ Co + c_ie + 0(e^ 



(A.IO) 



where 

Cl = 



1 

48 



log^ X — 6 log(x + 1) log^ X — 12Li2(— x) log x + Stt^ logx) 



- Grr^ log{x + 1 + 12Li3(-a;) - I2C3 



(A.ll) 



Co 



1 

960 



5 log'' X - 40 log(x + 1) log^ X + 120 log^(a; + 1) log^ x + lOvr^ log^ x 

- 1207r^ log(x + 1) logx + 480 log(x + l)Li2(-x) log x - 240 Li3(-a;) log x 

+ 4805i,2(-x)loga;- 240C3loga;+ 1207r2log^(a; + 1) -480 log(a; + l)Li3(-a;) 



+ 480Li4(-a;) - 48052,2(-x) + 480 log(a; + 1)C3 + vr^ 



(A.12) 
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2 log^ X + 10 log(x + 1) log^ X + 10 log^(x + 1) log^ X + 20Li2(-a;) log^ x 

- Stt^ log^ a; - 20 log^(x + 1) log^ x - SOCs log^ x + SOvr^ log^(a; + 1) logx 

- 1201og^(x + l)Li2(-x)loga; + 1201og(x + l)Li3(-x) logx - 120Li4(-x) logs 
+ 12052,2(-x) log(a;) - 240 log(x + l)5i,2(-x) log x - 2405i,3(-x) log x 

+ 120 log(x + 1)C3 logx + 47r^ log x - 207r2 log^(x + 1) 

- 87r^log(x + 1) + 1201og^(x + l)Li3(-x) 

- 40Li2(-a;)Li3(-a;) - 240 log(x + l)Li4(-x) + 240Li5(-x) + 240 log(x + l)52,2(-x) 
+ 40if2,3(-a;) + 120i73,2(-a;) + 40Li2(-x)5i,2(-a;) - 40/f2,i,2(-x) 



- 120if2,2,i(-2;) - 240C5 - 1201og2(x + 1)C3 - SOtt^C, 



(A.13) 



A. 5 The Half-Curtain Diagram 



We now consider the "half- curtain" diagram, whose contribution to the Wilson loop 



Vi 




is 



Figure 12: Diagram of the half- curtain topology. 



Whc(a;) = / dadpdri / rfr; 



Jti 



'2- 



— SaT2y^'' ( — STi — tp — UpTiY^'' 

{1 + x) ^ da da db dT2 



uj ' ' Jo Ji Ji Ji^ ' {ar^y^^ {1 - aby+^ ' 

(A.14) 

where we have changed variables in the second line to a = 1 + up/s, and b = l + uTi/t. 
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The evaluation of this diagram and of that with s ■H- 1 leads to 



Whc(x)+Whc(l/x) = -^f- 



u 



-2e 



X 



27recot(27re) - -x"^' 2Fi(2e, 2e; 1 + 2e; -1/x) 
^x2^2Fi(2e,2e;l + 2e; -x) 



+ 



(1 + x)'' + (1 + (^x-^ 2Fi(e, e; 1 + e; -1/x) 

+ 2Fi(e, e; 1 + e; — x) — 27re cot(7re) 
This can be expanded in e using 

2Fi(e,e;l + e;x) = 1 + e^Li2(x) - [Li3(x) - ^i2(x) 



(A.15) 



+ [Li4(x) - S22ix) + Si3(x)] + O(e^) , (A.16) 

The contribution of all diagrams of the half-curtain type is obtained by multiplying 
(lA.lSP by a factor of four. One obtains thus 



-s 



-2e 



+ i-t) 



-2e 



X 



+ do + d_ie + 0(e' 



(A.17) 



where 



di = 2ci 

do = - 

7 

2 
1 

96 



3co - ^ (Stt^ - log' x) (log' X + tt') 



(A.18) 
(A.19) 

c_i (A.20) 
- log^(x) + 6 log(x + 1) log^ X + 12Li2(-x) log^(x) - Svr' log^ x 



+ 67r' log(x + 1) log' X - 12Li3(-x) log' x - SOCs log' x - 427r'C3 
and Cj are the coefficients for the Y diagram, given in (IA.lip - (IA.13p . 



A. 6 The Cross Diagram 



We now consider the cross diagram, whose expression is given by 



JO JO 



1 /st\"'' . r\ r'",, i i 



^/ i-i/^K J-x Ji Ji (1 - aib2y+' (1 - a2biy- 



(A.21) 
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Figure 13: One of the cross diagrams. As before, the remaining three can be generated 
by cyclic permutations of the momentum labels. 



for which we find 

We,(a;) = 



1 /st 



-2e 



I - i 3F2(2e, 2e, 2e; 1 + 2e, 1 + 2e; -x)x^' 
r(-e + 1)2 3F2(2e, e, e; 1 + e, 1 + e; -x)x 



+ 



1 

+ 4 



r(-2e + ll 



2Fi(e, e; 1 + e; -x)x^ + 2F1 ^e, e; 1 + e; 



27re cotfeTr) 



X -H- - 

X 



e^7rcot(27re) {ip{2e) + 7) - vr^e^ cot(7re) 



(A.22) 



Notice the presence of the one-loop finite diagram squared. Expanding this and 
multiplying by a factor of two to account for all diagrams leads to the following 
result, 



\-2e 



fo + f-ie + 0{e' 



(A.23) 



26 



where 
/o 

/-I 



192 
1 



2 \2 



(A.24) 



2880 



3 log^ X + 30 log(x + 1) log^ X + 180Li2(-a;) log^ x 



X 



- 2071^ log^ X + eOvr^ log(x + 1) log^ x - 900Li3(-a;) log^ x - 540 Cs log^ 
+ 1807r2Li2(-x) log x + 2880Li4(-x) log x 

- 497rMogx + 307r^log(x + 1) - 4207r2Li3(-x) -4320Li5(-x) 

+ 432OC5 - 1020 vr^Cgj . (A.25) 



A. 7 The Hard Diagram 

The generic n-point hard diagram topology is depicted in Figure HI In the four-point 
case, one diagram is obtained from Figure H] by simply setting Q2 = Q3 = 0, and 
Qi = Pa- There are four such diagrams, obtained by cyclic rearrangements of the 
momenta. We have evaluated the hard diagrams using Mellin-Barnes, arriving at the 
following result: 



hard 



-2e 



X 



(A.26) 



where 



ho = — , 
48 ' 

hi = - ici 



1 
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ho = -2co - — (log^ x + n^) - -n^ (log^ x + n^) + ^^^^ 



vr 



(A.27) 
(A.28) 
(A.29) 
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200i/2,3(-a;) + 600i/3,2(-a;) - 200i72,i,2(-a;) - 600i/2,2,i(-a;) 



480 

+ 200C3Li2(-x) - 200Li2(-x)Li3(-x) - 200Li2(-x)Li3(x + 1) + 1320Li5(-x) 
+ 20Li2(-x) log^(x) + 60Li3(-x) log^(x) - 600Li2(-x) log^(x + 1) log a; 

- 600Li2(a; + 1) log^(a; + 1) logx + 100Li2(-x) log(-a;) log^(x + 1) 

+ 600Li3(-a;) log^(x + 1) + 20n^U2{-x) log{x) + 600Li3(-a;) log(x + 1) log(a;) 

- 600Li4(-a;) log(a;) + 1200Li4(x + 1) logx + 200Li2(-x)Li2(x + 1) log(2; + 1) 

- 1200Li4(-x) log(2; + 1) + 60052,2(-a;) logx + 120052,2(-a;) log(x + 1) 

- 240C3 log^ X - 6OOC3 log^ix + 1) - 6OOC3 log(x + 1) logx - 2 log^ x 

+ 5 log(a; + 1) log"^ x - 400 log(-x) log^(a; + 1) log x - lOOvr^ log^(a; + 1) 

- 407r^ log(x + 1) log^ X + ISOvr^ log^(x + 1) logo; + 50 log^(x + 1) log^ x 

- 100 log^(a; + 1) log^ x + Ttt^ logx - SSvr^ log(a; + 1) - IO2OC5 - 3207r=^C3 

(A.30) 

The analytical evaluation of this diagram up to 0{e^) was obtained in [1]. Our 
evaluation of the 0(e°) terms agrees precisely with that of [28] (and with [1] up to 
a constant term). The evaluation of the 0{e) term is new and has been performed 
numerically. We have then compared the entire Wilson loop expansion to the analytic 
expression for the amplitude remainder given in fl3.8p . finding the relation fl2.17p . 



B The One-Loop Wilson Loop Reloaded 



In this section we derive a new expression for the all-orders in e one-loop finite Wilson 
loop diagrams, and hence also a new expression for the all orders in e (finite part of the) 
2me box function. We also improve a previous expression for use in all kinematical 
regimes. 

A general one-loop Wilson loop diagram is given by the following integral 



r(l + e)e^^ X J-,(s,t,P^Q2) 
r(l + e)e^^ X J^' dr da- ^ 



{-[P^+a{s~P^)+T{t-Q'2+aTu)] ~ie}^+' ' 

(B.l) 

Here we have defined s = {p + Py, t = {p + QY and u = P'^ + Q'^ — s — t. The relation 
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between this Wilson loop diagram and the corresponding 2me box function is [3] 

y \ ,r(l + e)r2(l-e 



finite part of 



rfl-2e) 



(B.2) 



Notice that we have included an infinitesimal negative imaginary part —ie in the de- 
nominator which dictates the analytic properties of the integral. This has the opposite 
sign to the one expected from a propagator term in a Wilson loop in configuration 
space. On the other hand it has the correct sign for the present application, namely 
for the duality with amplitudes One simple way to deal with this is simply to 
add an identical positive imaginary part to all kinematical invariants 



s s + ie, t t + ie, 
We will assume this in the following. 



^P^ + le, -^Q^ + ie . 



(B.3) 



Changing variables to cr' = u — (P^ — t)/u and r' = r — (P^ — s)/u and then 
dropping the primes, this becomes 



T,{s,t,P\Q^^ 



u/2 



da I dr — . 

2 , 



1 
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da 



2e yt-p£ a [-(a-i + crrw)]' 



T = {Q^-t)/u 



^{s~P2)/u 



(B.4) 



(B.5) 



where a = u/{P'^Q'^ — st). Note that the analytic continuation of a implied by ( IB. 31) 
is a ^ a — ie. Now we split the integration into two parts. 



+ 



(B.6) 



and rescale the integration variable so that it runs between and 1 in each case. It is 
important to split the integral in this way since there is a singularity at a = which 
one must be very careful when integrating over. We obtain in this way 



2e 



1 



a 111 
1 



[l-aP^)aY [l-il-aQ^)a] 
1 



[l-{l-as)aY [l-il-at)aYj ' ^^'^^ 

where a = u/ {P^Q^ — st). Now each of the four terms by itself is divergent (even for 
e 7^ 0), only the sum gives a finite integral. A straightforward way of regulating each 
term individually is to simply subtract 1/a from each term in the integrand, thus 
removing the divergence at o" = 0. We thus have 



^.(s,t,p^Q' 



/(l-aP2) + /(l-aQ2)_/(i_as)-/(l-at) , (B 
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where 

f{x) = - / da- = - / da- . (B.9) 



^ Jo ^ Jo or 

The problem becomes that of finding the integral f{x). It has two equivalent forms, 
both given in terms of hypergeometric functions. The first form is given by 

oo 

fix) = xx3F2(l,l,l + e;2,2;x) = ^e"^i„+i(x) . (B.IO) 

n=l 

Notice the very simple expansion in terms of Nielsen polylogarithms. The second 
form is 

fix) = [(-x)~^2Fi(e,e;l + e;l/a;) + eloga;] + constant , (B.ll) 

where the constant is there to make /(O) = 0, and is not important since it will cancel 
in 

We thus arrive at two different forms for the Wilson loop diagram. The first form 

is 

T,is,t,P^Q') 



-a) 



[1 - aP%F2il, 1, 1 + e; 2, 2; 1 - aP^) + (1 - aQ%F2il, 1, 1 + e; 2, 2; 1 - aQ^) 



- (1 - as)3F2(l, 1, 1 + e; 2, 2; 1 - as) - (1 - at)3F2(l, 1, 1 + e; 2, 2; 1 - at) , 

(B.12) 

and it is manifestly finite. Furthermore, since 3F2(1, 1, 1 + e; 2, 2; x) = Li2(x)/x, this 
form directly leads to the expression derived in [IHIHT] for the finite 2me box function, 

J-,=o(s,t,P',g') = ^ [Li2(l-ap2) + Li2(l-aQ2)-Li2(l-as)-Li2(l-at)] . 

(B.13) 

We also notice that the simple expansion of (IB.IOI) gives a correspondingly simple 
expansion for the Wilson loop diagram in terms of Nielsen polylogarithms. 

The more familiar looking second form for the two-mass easy box function is (see 
(A. 13) of HH]) 

T,is,t,P^Q') =--^x 



a 



1-aPy " ' ' ' " \l-aQ 



2Fi(e, e; 1 + e; 1/(1 - aP^)) + — 2Fi(e, e; 1 + e; 1/(1 - aQ^)) 



2Fi(e, e; 1 + e; 1/(1 - as)) - ( ) 2Fi(e, e; 1 + e; 1/(1 - at)) 



1 — as/ \l — at^ 

+ ei-a)' log(l - aP2) + log(l - aQ^) - log(l - as) - log(l -at)] . (B.14) 
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This second form was derived in [SlHH] except for the last hne which is an additional 
correction term needed to obtain the correct analytic continuation in all regimes. The 
identity 

(l-an(l-aO=) _ ^ ^^ ^^j 



^1 - as){l - at) 



implies that if all the arguments of the logs are positive then this additional term 
vanishes, but for example if we have 1 — aP^, 1 — aQ"^ > and 1 — as, 1 — at < 
then the additional term gives (taking care of the appropriate analytic continuation in 
(IB. 311 ) sgn(a)27ri(— a)'^/e. This becomes important when considering this expression 
at four and five points in the Euclidean regime. 

For applications in this paper we are interested in taking either one massive leg 
massless (for the five-point case) or both massive legs massless (for the four-point 
case). Using the first expression for the finite Wilson loop diagram in terms of 3F2 
functions and using that 3F2(1, 1, 1 + e; 2, 2; 1) = ""^(^"^^"^ = — (where ^p{x) is the 
digamma function, 7 is Euler's constant and if„ is the harmonic number of n), we 
obtain the four- and five-point one-loop Wilson loop expressions of (14. ip and (14.21) . 

For completeness we also consider the limit with = = using the second 
expression for the finite diagram (lB.14p . since this and similar expressions have been 



used throughout appendix A. When P"^ = = 0, we have a = 1/ s + 1/t, 1 — as 
—s/t and 1 — at = —t/s and using 



2Fi(e,e;l + e;l) =e7rcsc(e7r) , (B.16) 



we get 



J-,(s,t,0,0) = -^x 



^ (S)' (j)' ^ + -'/^) -©'(?)' ^ + 

+ 2(a)'e7rcsc(e7r) +e(-a)'(|- log(l- as) -log(l- at)) . (B.17) 

We wish to know this in the Euclidean regime in which s,t,a < 0. The first line is 
then manifestly real, whereas the second gives 

27re(-a)^e-^'^ csc(e7r) + 27rei(-a)' = 27re(-a)' cot(e7r) . (B.18) 

This is the form used for the one-loop Wilson loop throughout the paper, for example 
in dH]) and (IX5|) . 
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